Quantum Rings with Rashba spin orbit coupling: a path integral approach 
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We employ a path integral real time approach to compute the DC conductance and spin polariza- 
tion for electrons transported across a ballistic Quantum Ring with Rashba spin-orbit interaction. 
We use a piecewise semiclassical approximation for the particle orbital motion and solve the spin 
dynamics exactly, by accounting for both Zeeman coupling and spin-orbit interaction at the same 
time. Within our approach, we are able to study how the interplay between Berry phase, Ahronov 
Casher phase, Zeeman interaction and weak localization corrections influences the quantum inter- 
ference in the conductance within a wide range of externally applied fields. Our results are helpful 
in interpreting recent measurements on interferometric rings. 

PACS numbers: 03.65.Vf,72.10.-d,73.23.-b,71.70.Ej 



I. INTRODUCTION 

In classical physics, a charged particle moving in an 
external magnetic field B feels a (Lorentz) force only in 
the regions in which B is different from zero. Instead, 
in 1959 Ahronov and Bohm (AB)^ showed that the wave 
packet representing the state of a quantum particle can 
be influenced by an external vector potential A, even if 
the corresponding magnetic field is zero, provided that 
the particle is moving in a space with a nontrivial topol- 
ogy (holes). To be more specific, the wavefunction of the 
charged particle may acquire a nonzero phase, when un- 
dergoing a closed path in a space threaded by an external 
magnetic flux. 

Nowadays, mesoscopic quantum rings (QR's) allow to 
have direct access to the phase of the electron wavefunc- 
tion, since their size is smaller than the distance over 
which the phase randomizes, as a consequence of scatter- 
ing against impurities and interactions. The total rate 
of particles coming from the two arms of the ring and 
interfering at the exit contact can be directly probed by 
measuring the QR conductance. Indeed, interference ef- 
fects have been observed in metal QRs many years ago^. 
In addition, since electrons are spinful particles, the spin 
part of the wavefunction is influenced by the magnetic 
field via the Zeeman term in the Hamiltonian. Moreover, 
a more subtle effect arises when there is a magnetic field 
non-orthogonal to the plane of the orbiting particle be- 
cause, as a consequence of the orbital motion, its spin dy- 
namics is instantaneously governed by a time-dependent 
Hamiltonian'^. This time dependence ends up in an ex- 
tra phase acquired by the particle wavefunction which is 
named after Berrj^^S., who put in foreground its topolog- 
ical properties when the orbits are closed. 

Recently, semiconductor technology allows to grow 
samples (for instance made by InAs or InGaAs) with size- 
able spin orbit interaction (SOI). Rashba^ first pointed 
out that the SOI strength can be controlled by means 
of voltage gates. This feature has been recently found 



experimentallyiiS,. Because of the SOI, an electric field 
E orthogonal to the orbit provides a momentum depen- 
dent effective magnetic field felt by the electron spin in 
orbital motion. Such a field can be equally well described 
by a vector potential that adds an extra phase to the 
wavefunction^. This phase was spelled out by Ahronov 
and Casher-iS, as a "dual" AB effect, with charge and 
spin interchanged (together with B and E). 

During the last years, the effects of SOI on the AB 
oscillations have been observed in semiconductor based 
QR by several groupsii»i2,il^Jil. It has been recognized 
that, in the presence of both B and E orthogonal to the 
orbit plane, the effective, momentum dependent total B 
field is tilted w.r. to the vertical direction. The resulting 
Berry phase influences the interference pattern. Such 
a result is of the utmost interest because the Rashba- 
SOI (RSOI) turns out to be a tool to tune the Berry 
phase and, ultimately, the conductance, as well as the 
spin polarization of the outgoing electrons. In RefJ^, it 
is clearly shown that the AB interference fringes can be 
modified by tuning an external electrostatic potential. 

In a recent publicatioi>i^ we have included all these 
phenomena affecting the interference of an electron bal- 
listically transported in a ring, also accounting for some 
dephasing at the contacts. We have also shown that the 
AB peak in the Fourier transform of the magnetocon- 
ductance displays satellite peaks due to the RSOI. It has 
been suggested that satellite peaks observed in recent 
experimentsii have the same origin. In the present work 
we review the theory by presenting the full calculation. 
In addition, we include the semiclassical paths leading to 
weak localization corrections^^, and we discuss the rota- 
tion of the spin polarization, during the transport across 
the ring. 

During the last years, several theoretical techniques 
have been employed to study quantum transport in a QR. 
In Ref.s?^?-, an imaginary time path integral approach- 
is developed to study the conductance of a strictly one 
dimensional (ID) QR and its conductance fluctuations 
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in the diffusive limit. In Ref.^^ a real time path inte- 
gral approach is applied in the limit of negligible Zeeman 
splitting. Several papers have discussed the conductance 
properties and the spin selective transport of QR's in 
the strictly ID ballistic limit, by means of a spin depen- 
dent scattering matrix approacb22iS^i2iii2^i2^. In the ab- 
sence of the magnetic flux, the conductance shows quasi- 
periodic oscillations in the SOI stength, which can be 
modified by switching the magnetic field on. Numerical 
calculations^li^i^ have shown that in the 2D case there 
are only quantitative modifications of the ID results that 
do not qualitatively affect the physics. 

In this paper we extend the real time path integral ap- 
proach previously developed in Ref.— to study the con- 
ductance and the spin transport properties of a ballistic 
quantum ring in the presence of both RSOI and of an 
external magnetic flux orthogonal to the ring plane. We 
use a "piecewise" saddle point approximation for the or- 
bital motion, keeping the full quantum dynamics of the 
spin. This approach allows us to take into account, in a 
nonperturbative way, both the RSOI and AB phase and 
to include also the Zeeman spin splitting. 

Our numerical approach evaluates all paths contribut- 
ing to the quantum propagator. The scattering at the 
leads can be forward or backward, according to the prob- 
ability amplitudes given by the S-matrix. Weak localiza- 
tion corrections can be easily extracted from our result. 
We also allow for some diffusiveness at the contacts by 
adding a random phase factor in the motion. 

The DC conductance is derived from the Landauer 
formula^o 5 = /TiY^aa' \A{a;a'\E)\^ , where A{a]G'\E) 
is the probability amplitude for an electron entering the 
ring with energy E and spin polarization a' to exit with 
spin polarization a. We also report the change in the spin 
polarization the electron transported across the ring. 

The structure of the paper is as follows: 

• In Section II we introduce the Feynman propagator 
for a spinful electron injected at the Fermi energy 
in the ring. 

• In Section III we discuss the topology of the allowed 
paths and the scattering of the electron at the leads. 

• In Section IV we represent our path integral in the 
coherent spin basis^i^^ and derive the saddle point 
equations of motion, whose classical counterpart is 
described in detail in Appendix A. This allows us 
to justify the choice of a piecewise semiclassical ap- 
proximation for the orbital motion of the electron 
in the ring. 

• In Section V we present the details of the calcula- 
tion by rewriting the path integral as a collection 
of single arm propagators. These are the building 
blocks to be calculated in the next section. 

• In Section VI we analyze how the orbital motion 
affects the full quantum dynamics of the electron 
spin for each arm of the ring and chirality. The 



spin propagator is derived in Appendix B, in the 
basis corresponding to the rotating reference frame 
in the spin space. 

• In Section VII we discuss the dependence of the 
conductance on the external fields and on the over- 
all transmission across the ring. 

• In Section VIII we focus on the spin polarization of 
the outcoming electron. 

• Section IX includes a short summary and our con- 
clusions. 

II. THE TRANSMISSION AMPLITUDE 

Our model Hamiltonian will be the two-dimensional 
Hamiltonian for a particle with spin-1/2 5, in an or- 
thogonal magnetic field, with spin-orbit coupling to an 
orthogonal electric field (Rashba coupling). It is given 
by 

H[p, fj] = ^ (p+ -Ao) ' - We 5z + Hso (1) 
2m \ c / 




where a is the spin orbit coupling constant, in units 
eV K, S — ha/2 {<Tx,(Jy^az are the Pauli matrices), 
^o(^) = -f(— y,^;, 0) is the vector potential generating 
the uniform field B, normal to the ring surface, taken in 
the symmetric gauge, Uc — geB/2mc is the cyclotron fre- 
quency. In real nanostructures based on InAs or InGaAs 
the g factor can strongly deviate from the value of two. 
However the result we present here, are fully general as 
they depend on the ratio a/hu>cR which can be tuned by 
acting on a. Since we will assume only a single channel 
to be available for electron propagation across the ring, 
we will picture the single channel ring as a 1 — d circle 
of radius i?, connected to two leads. Accordingly, the 
position of the particle within the ring is parametrized 
by the angle ip and the vector potential has just the az- 
imutal component = (j)/2TTR, where 4> is the magnetic 
flux threading the ring. 

In order to study the conduction properties of the ring, 
one needs the propagation amplitude for an electron en- 
tering the ring with spin polarization /io to exit with spin 
polarization /iy, at energy Eq. This is given by 

f°^' dtf - Eq'/ 

A{nf;fio\En) ^ — " (^/,M/,t/|n),A'o,to), (2) 
Jo To 

where (f/, i/|ro, /iq, io) is the amplitude for a par- 
ticle entering the ring at the point fo and at the time to 
with spin polarization /io to exit at the point ff at the 
time tf with spin polarization fXf. In our tensor product 
notation, we define \r,fi) = \r) . tq — mR'^/{2h) 

is the time scale for the quantum motion. 
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In order to compute {ff,^f,tf\fo,iio,to) , we resort 
to a path integral representation for the orbital part of 
the amplitude. Since we parametrize the orbital motion 
of the particle in terms of the angle ip, we provide the 
pertinent Lagrangian, Corb, as a function of (p, if. It is 
given by 

(3) 



The last two contributions to Eq.Q are a constant, com- 
ing from the spin-orbit term, and the Arthurs^ term, 
which is required when a path integration is performed 
in cylindrical coordinates. Since both contributions are 
constant, they can be lumped into the incoming energy 
Ea and therefore they will be omitted henceforth. By 
taking into account the spin degree of freedom, as well, 
we represent the propagation amplitude as 



I 



where q — (t>/4>o, ipo being the flux quantum hc/e. 



Uspin{tf,h) = Texp 



dt' Hspinit') 



(5) 



is the full spin propagator and the spin Hamiltonian 
is given by 



^ H spill (i) 



^<^e*¥'(*) -f 



(6) 



with 7 = 2aTQ/{hR)^. In Sections IV and V we show 
that the amplitude of Eq.(|3]) can be approximated by 
choosing a piecewise semiclassical orbital motion for the 
particle in each arm of the ring, while keeping the full 
quantum dynamics of the spin. In particular, we will see 
that, within the physically relevant range of parameters, 
the orbital motion can be approximated as a uniform ro- 
tation (with constant angular velocity), which makes the 
spin dynamics to be the one of a spin-1/2 in an effective, 
rotating, external magnetic field. Yet, in order to explain 
how we deal with quantum backscattering at the contacts 
between ring and arms and corresponding dephasing ef- 
fects, we will introduce our formalism in the next section, 
by discussing a simplified version of our problem: a spin- 
less electron propagating across a mesoscopic ring. 



III. FEYMAN'S PATHS FOR A SPINLESS 
PARTICLE TRANSMITTED ACROSS A RING 

In this section we introduce our formalism by comput- 
ing the transmission amplitude for a spinless electron of 
mass m and charge — e, traveling across the ring in an 
orthogonal magnetic field. For a realistic device, at each 
lead one has to take into account three possible scatter- 
ing processes, consistently with the conservation of the 
total current. This is described in terms of a unitary 
5— matrix that, in the symmetric case in which the two 



arms are symmetric, is given by 

/ -i(l + f) 1(1 -f) 
i(l-r-) -i(l + f) 



S 



F2) \ 

r2\ 



1 



(7) 



The numerical labeling of the S-matrix elements referring 
to the three terminals of each contact fork, are explained 
in Fig.([Tl la). Assuming, for simplicity, that the scatter- 
ing matrix is the same for both leads, Eq.([7|) will hold 
both at the left-hand lead, and at the right-hand lead of 
the ring. 

In particular, 53^3 = f is the reflection amplitude for 
a wave coming from the left lead, iSi(2),i(2) = — ^(1 + ^) 
the reflection amplitude for a wave incoming from the up- 
per/lower arm, 5i(2),2(i) = ^(1 ~ ^) is the transmission 
amplitude from the upper (lower) to the lower (upper) 

arm and 5i(2),3 = '^3 1(2) = y^l(l — f^) is the transmis- 
sion amplitude from the upper/lower arm to outside of 
the ring. 

In Fig. 11]) we show the simplest paths of the electrons 
in the ring including forward scattering at the contacts, 
only. The contribution to the total amplitude coming 
from such paths, in which the electron enters the ring at 
an angle (po at time to and exits at ipf at time t/ , is given 

by 

A{p^f,tf;po,to)= E / P¥.(r) e-^-^Mt)]/^ 

= g-.g(¥>/-v.+2™) / p^(i)e-*'^-/'*o '^*'^''(*) ,(8) 

where we have summed over paths in which the electron 
winds n-\-l/2 times in the ring, before exiting it. Positive 
(negative) n values imply clockwise (counterclockwise) 
propagation along the ring. 
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This propagator can be evaluated exactly^^. However 
we report here just the saddle point evaluation, for com- 
parison with the spinful case. Minimizing the action gives 
the classical equation of motion (together with the per- 
tinent boundary conditions for a path winding n -\- 1/2 
times): 



Let us assume that the particle is injected in the ring at 
(/7o = and comes out a,t ipf = tt in a time T = tf. Eq. ([H]) 
gives: 

I +00 ' 

V nitf ^ 

' n— — oo 

where the prime in the sum takes into account the fact 
that one does not sum over rt = 0, and the square root at 
the prefactor accounts for the gaussian fluctuations. Of 
course, this propagator is periodic in q of period q = \ 
up to a minus sign. 

More involuted paths arise if one takes into account 
backscattering processes in which the electron can get 
backscattered within the same ring's arm from which it 
is coming. For instance, the paths (2/) and (2/i), as 
well as (2(7) and (2i) in Fig.©, include looping in op- 
posite directions around the ring's hole. Interference be- 
tween clockwise and counterclockwise windings leads to 
weak localization corrections. We denote these corre- 
sponding paths (including also (2c) and (2c?) ) as "re- 
versed paths" (i?P). In our approach, all order paths 
are numerically generated up to the convergency and 
the S'— matrix of Eq.(l7]) is implemented in the numer- 
ical algorithm. In conclusion, we have established 




FIG. 1: (color online) First and second order paths included in 
the calculation of the transmission amplitude across the ring, 
from left to right, including forward scattering only. Numbers 
1, 2, 3(1', 2', 3') in Fig. la refer to the labeling of the terminals 
in Eq.©) 




FIG. 2: Second order paths of the transmission amplitude 
from left to right including backscattering at the leads. Paths 
(2/) and (2/i), as well as (2(?) and (2j) contribute to the weak 
localization corrections. 



that the paths contributing to the transmission ampli- 
tude can be built up by adding four types of elementary 
paths: forward orbiting in the upper arm of the ring 
(</? G (0,7r)), u<_ backward orbiting in the upper arm of 
the ring [ip G (tt, 0)), forward orbiting in the lower 
arm of the ring ((/? £ (27r,7r)), backward orbiting in 
the lower arm of the ring {ip e (7r;27r)). In Section V, 
we will generalize such an approach, while in the next 
Section we discuss the Feynman path-integral represen- 
tation of the propagation amplitude for a spinful electron 
propagating along one of these elementary paths. 



IV. QUANTUM AMPLITUDE AND 
SEMICLASSICAL ORBITAL MOTION FOR A 
SPINFUL ELECTRON IN THE RING 

In this section, we construct the Feynman represen- 
tation in the basis of the coeherent spin states, for the 
propagation amplitude of an electron moving along one 
of the arms of the ring with a given chirality. To discuss 
the saddle-point approximation we need the equations of 
motion coming from the condition that the action is sta- 
tionary. The coherent spin state basis provides a straight- 
forward route to perform a semiclassical approximation 
involving both orbital, and spin degrees of freedom, at 
the same time. In particular, we will show under which 
conditions the classical orbital motion = cnst, can be 
retained, as in the spinless case. Accordingly, the spin 



(p{t) = ; Lp(ti) = Lpi , ip(tt) = ipf + 27rn 



(9) 
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dynamics will be that of a "quantum magnetic moment" 
in a time-dependent external magnetic field. 

Let il denote the orientation of the spin S and be 
the coherent state such that: 



for the three components of the spin vector, respectively. 
The full propagator in the coherent spin state represen- 
tation is given by: 



{n\s\n) = s[n] = ns 



sin Q cos $ 
sin 8 sin $ 
cos Q 



(11) 



where: 



S. 



spin 



ve / I?$ ^-^s.,.40,e;<fM^,^] ^ (12) 

e(o)=eo J*(o)=*o 



(l-cosO) 



(13) 



with: 



6, 6; $, ^\(p, if = cos 6 — 793sin6 cos((/7 — 4') 



(14) 



r 



The Lagrangian of Eq. (|T4[) corresponds to the coherent 
spin-state representation of the classical Lagrangian de- 
rived from Eq.([T]) 



In order to extract the relevant physics from the above 
equations, we multiply Eg. p^ by ip^ and, by use of 
Eqs. (|17ll8p . we rewrite Eq.ldl]) as: 



-^b, ^, s] = 



• A + 2ma 



X S 



- LUcSz .(15) 



In Eq. (fT5)) we have introduced the constraint that the or- 
bital electron motion takes place along a one-dimensional 
circle by parametrizing the trajectories with the angle tp 
as X = Rcosip; y — Rsimp. The additional Berry phase 
termi2£ in Eg. p^ arises from the fact that different spin 
coherent states are not orthogonal to each other since, 
to leading order in e, the scalar product between spin- 
coherent states at times tj,tj + e, \^{tj)) and \^{tj + e)) 
is given by 



{n{tj + e)\n{tj)) w exp 



-[1 - cose{tj)]^{tj)e 



(16) 



Let us look for the trajectories in orbital and spin space 
which make the action stationary. This requires solv- 
ing the Eulero-Lagrange equations for the Lagrangian 
C[(p, (p; 8, 8; <&, <i>] appearing in Eq. fT^ . which are given 

by 



sin 8 -|- 793 sin x sin 8 



sin8[-<l> 



0, (17) 
27(^9 cos 8 cos X — 0, (18) 



(/3 - — [8 cos 8 cos X + <i> sin 8 sin x] = 0, (19) 
with X — ^ ^ ^■ 



ifcp = — 8 sin 8 . 

A straightforward time integration gives: 



— cos 8 — cnst 
4 



(20) 



(21) 



which states that the total energy is conserved. Ac- 
cording to Eq. pTjl . the particle energy only includes the 
orbital kinetic term and the Zeeman term. This is obvi- 
ous, since the force associated to the spin orbit coupling, 
being gyroscopic, does no work. A change in the preces- 
sion angle 8 implies acceleration in the orbital motion. 
According to the spin Hamiltonian of Eq.® the RSOI is 
responsible for flipping of the spin, while the Zeeman cou- 
pling tends to stabilize the spin direction. Hence, there 
are two physically relevant limits, in which the orbital 
motion fully decouples from the spin dynamics, accord- 
ing to the inequality lOc/2 < (>)<^, as we are going to 
show next. Both limits involve an orbital motion with 
constant velocity ip and a spin orientation given by the 
angles 8 ~ constant , — $ = {mod. tt). 

a) Vanishing Zeeman coupling: This case has already 
been considered and it has been shown that it can be ex- 
actly solved analytically^!. In this case Eq. (PO)) allows for 
the semiclassical solution ip = cnst and quantum fluctua- 
tions of the spin do not interfere with the orbital motion. 
The spin is tilted by the constant angle tan 8 = 27 and 
precesses with a constant frequency $ = </?. 
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b) Large Zeeman coupling: flipping of the spin 
and quantum fluctuations of the spin are strongly dis- 
favoured, again corresponding to the classical saddle 
point configuration O = constant , — $ = 0, tt with 



n/tf — If — ±<I> = LUc 



1 



1 — 27 ctanO 



(22) 



The spin precesses around the z-axis, with the same an- 
gular velocity as the orbital motion. When 7 ^ 0, 8 can 
be vanishingly small, with ^ — ujc- 

It may appear that the simple spin precession provided 
by Eq.® with 6 = cnst is a solution of Eq.s (|17ll8ll9p 
for any ratio uJc/'f- Were this the case, there would be no 
need to invoke the limitations of case b) as stated above. 
However, a careful analysis of the stability of this saddle 
point solution shows that the simple spin precession is a 
minimum of the action only when the Zeeman coupling 
is strong. In particular, if Uc/f >> 1, an additional 
condition lOcToIi^ > 1 has to be satisfied. Otherwise 
the frequency of the fluctuations around the saddle point 
solution does not keep real and the analysis of small os- 
cillations in the parameter space around the saddle point 
solution breaks down. 

This shows that a classical orbital motion with con- 
stant velocity is compatible with both limits of large and 
small ratios of the Zeeman coupling to the RSOI strength. 
However the two limiting cases do not belong to the same 
saddle point. In particular there will be a crossover region 
connecting the two limits in which quantum fluctuations 
of the spin may induce changes in the orbital velocity of 
the particle. 

In the rest of the paper we will choose tp — const 
piecewise and parametrize the quantum dynamics of the 
spin with the value of the velocity obtained by the sta- 
tionary phase condition discussed in the next section. 
As discussed above, our approximation cannot repro- 
duce faithfully the intermediate region of parameters 
ranging the two limits of large and small Zeeman cou- 
pling ujc/'y'P < (>)!• We checked numerically the reli- 
ability of our approximation, by numerically integrating 
Eas. (|17ll8ll9p and have have found that it is satisfactory 
in most of the parameter range. 

By putting ip = cnst into Eqs. (|17ll8p . they be- 
come, as we show in detail in Appendix A, the clas- 
sical equations of motion for a magnetic moment in 
a time dependent magnetic fleld B = {B+,B-,Bz) — 
[jif e'"^, 7(^ e~^'^, —ujc/i) ■ This can be easily understood 
from the fact that minimizing the action in Eq. p^ di- 
rectly w.r. to the spin components provides: 



S X 



'dt 



5H 



(23) 



which has to be solved together with the constraint of 
constant modulus: S ■ dS / dt — ( see Eq. (jA6|) . 

Among the other possible saddle-points, a solution of 
the motion equations can be found with the particle 
trapped within the ring arm (turning points of ip{t) are at 



((9 = 0, tt) if 7 is large enough. This solution doesn't seem 
to be practical as it requires fine tuning of the external 
parameters with transfer of energy from the spin motion 
to the orbital motion and viceversa. We have not inves- 
tigated it in detail, but we expect that it could provide 
resonant tunneling across the ring. 



V. SADDLE POINT APPROXIMATION AND 
LOOPING IN THE RING 

In this Section we show how we implement the piece- 
wise saddle point solution for the orbital motion, ip = 
cnst, to study the coherent propagation of the electron 
inside the ring. In the following, we will denote by 
"loop" and "looping trajectory" both trajectories that 
wind around the ring (closed), and paths in which the 
particle moves forth and back in one of the ring arms 
(open) (see Fig.®). Of course, the amplitudes differs 
very much in these two types of looping . Indeed, the 
net spin rotation is different between the two paths and, 
also, the Ahronov-Bohm {AB) phase is absent in the lat- 
ter ones . 



n=l 









i 




















An=3 



FIG. 3: (color online) Cayley tree describing the way in which 
higher order paths are built in the numerical code. Full lines 
correspond to propagation in the upper arm { u ) or lower 
arm ( d) in the forward (^) or backward ( <— ) direction. 
The nodes represent the leads of the ring. The exit nodes 
are marked by a broken line at each order n. The heavy line 
correspond to the path reported in Eq. |25|| . 



In general, at order n, we will have 2(^l"l~^^ trajecto- 
ries of a particle entering the ring at if — and exiting 
at (f — TT and each of them will include 2\n\ — l{n = 
±1, ±2, ...) elementary paths (or "stretches") of the type 
, , d^ , d^ , as classified at the end of Section II. 
Looping at order \n\ involves 2|rt| scattering processes 
at the contacts. In the case of spinful electrons the S*- 
matrix is doubled (6 x 6) w.r. to the one given in EqjT] 
Of course, the S'-matrix at the contacts is sample depen- 
dent and any special choice is arbitrary. In the following 
we neglect possible asymmetries in the up-down channel 
as well as accidental spin hipping in traversing the con- 
tacts. Hence, we make the simplifying assumption that 
the S'-matrix is block diagonal of the form given in Eq[7| 
for both spins. 

Each time the trajectory impinges at a contact with- 
out leaving the ring, the S— matrix of Eq.© alledges 
for either forward, or backward scattering. Let us de- 
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note with u^(^^){ti,tj) the forward (backward) propaga- 
tion ampHtude in the upper arm from time tj to time ti 
{d^(^^^{ti,tj) for the lower arm). The u{d) amphtudes 
include the Ahronov-Bohm phase and the spin evolution, 
but not the dynamical phase, which is factored out (see 



below). As shown in section II, to first order (|n| = 1) 
there are only two possible paths (see Fig. ([1] la, lb)), the 
propagation amplitude is, then, the sum of the the cor- 
responding two amplitudes: 



dtf .,fov , ^ , 
To 



(24) 



The amplitudes have to be summed all together, order by order. The key observation appearing in the symbolical 
writing of Eq. (j25p is that the dynamical phase, at a given order, does not depend on the chirality of the motion and 
can be factored out. On the contrary, the Ahronov-Bohm phase depends on the chirality, and the spin evolution 
depends on both the chirality and on the modulus of the propagation velocity. Beyond the first order, there is a net 
increase in the number and type of the paths to be summed together. In Fig. ([3]), we pictorially sketch all the possible 
paths by mean of a Cayley tree. Each node represents a lead of the ring and, according to the scattering matrix of 
Eq.©, the electron can be backreflected into the ring's arm it is coming from with an amplitude Sa, or, with an 
amplitude Sij, it can be either transmitted to the other arm, or outside of the ring. In the tree, the transmission out 
of the ring is possible at all the nodes crossing the dashed lines. Each dashed line is labeled by the order n of the 
interference in the ring. As an example, we explicitely write down one of the possible second order paths (the bold 
red hne marked by V in Fig. ((21), which corresponds to Fig.([2l2h)): 



^^{^J'f^■,^J'Q\Eo) 



^dtf^ 
tq 



tf rt2 

dt2 / dti{fif,tf\S3>2d- 

3 '^in 



.(t/,t2)522rf^(i2,ti)52iw^(ii,io)5i3|/^o,Ve"''^^*o'^*'^'^*'(25) 



At the saddle point with uniform velocity, (p = 2tt {2\n\ — l)/tf, we classify the collection of sequences belonging to 
a certain n as {C„}, and denote the superposition of amplitudes (e.g., the ones in the big parenthesis of Ea. ([M)) to 
first order ) as [fXf , ^o;q, ip\Cn\- In this way, we can rewrite the full amplitude of Eq.© as: 



A{fif;no\Eo) = 



dtf 



-iEotf /h 



+00 ' 

— y y 

n— — oo {Cti} 



T mo; 9, ^ - 2^ (2|n| ~ 1)A/|C„] ^(2|"l-i)^-o/t/ . (35) 



r 



The series in Eq. (j26|) is uniformly convergent. Thus, we 
may swap the integral with the sum, and integrate term 
by term. The integral contributing to order n is given by 



dtt 



'^0 g-J-Eot//?i-«7r^(2|n|-l)Vo/t/ 



= J- rf^g-.-^— ^(2H-i)V-^ (27) 

V TTl Jo 

with e — EqTq/H. We compute it approximately 
within stationary phase contribution. Since the phase 
of the exponent of the integrand is stationary at in — 
£(-1/2) 7^ (2\n\ — 1)to, by inserting this value in the phase 



and integrating out the gaussian fluctuations, we readily 
get: 



/ W e-»\/^27r(2|n|-l) 



W2e 



d{6x) e-'" 

-iv/i27r(2|n|-l)2 _ (28) 



(— ie^^/^ is the usual factor appearing in the one- 
dimensional free particle Green's function in real space 
and energy). This approximation, when plugged into 
Eq. ([26]) . provides the final result: 



+00 



A(M/;Moii?o) = ^ y (y :f 



fif,tn,fio;q,'P = 27r (2|n| - l)/t,; 

I 



-iVe27r(2|n|-l) 



(29) 



The enumeration of the trajectory configurations belong- ing to the collection C„, to order n, is numerically per- 
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formed order by order. 

In the next Section, we discuss the elementary prop- 
agators for each of the four stretches, u^,u^,d^,d^, 
as defined in Section III. This allows us to construct the 
functional for each incoming and outgoing spin polar- 
ization. 



VI. QUANTUM SPIN DYNAMICS OF THE 
ELECTRON PROPAGATING IN THE RING 

In this Section, we provide the explicit formula for 
the spin contribution to the total propagation amplitude, 
given by 



Uspin{t,t') ^ Texp 



(30) 



As discussed in detail in appendix A, within the saddle 
point approximation, ifspin(i) corresponds to the Hamil- 
tonian of a spin-1/2 in a time-dependent external mag- 
netic field. It can be written as (from now on, we will 
denote by ujo the frequency of the orbital motion, that 
is, the stationary phase value of (p) 



Thus, the matrix diagonalizing Hspin(t) at time t is 



cos I -sin|e~'"°* 
sin f e*"°* cos f 



(34) 



The matrix B{t) encodes the adiabatic contribution to 
the evolution of |^'(t)). Therefore, in order to write down 
the Schrodinger equation with Hamiltonian H, 

d 



i--H{t)\\^>{t))=0 



in the adiabatic basis, one has to strip off from the state 
\^{t)) its adiabatic evolution, operating with &{t), so to 
get: 



^^^-B^{t)H{t)B{t)+B{t)h^^B{t) 



B\m{t)) = Q. 

(35) 

Eq. (j35p may be rewritten in a 2x2 matrix formalism. 



Let 



be the components of |^'(t)) in the adiabatic 



basis. The corresponding system of differential equations 
reads 



Hsp\n{t) 

with 
rcosi? 



rcosd rsnv-de 
■sini9e*'^°* —rcosd 



UJc 

Y 



rsini? 



tan I? 



27a;o 



(31) 



(32) 



dt 



where we have defined 



Ha = 



a Ait) 



■WoSin^l icjosini^e"''^"* 



^ujo sin i9e 



(36) 



(37) 



Including only the AB phase implies d — 0, while in- 
cluding only RSOI implies "d ^ tt/2. 

It is useful to solve for the spin dynamics in the rep- 
resentation of the instantaneous eigenstates of Hspin{t). 
At fixed i, its eigenvalues are given by ±e = ±r, while 
the corresponding eigenvectors take the form: 



sin I e*'^°* 



\-,t) = 



-sinf e^*'^"* 
cos J 



The extra term appearing on the diagonal w.r.to the 
hamiltonian of Eq. pT]) is just the Berry phase: 



{+,t\i 



dt' 



,t) = -{-,tVj^\-,t) 



uj„ sm 



(38) 



Eg. ([55]) is solved in Appendix B and the full propagator 
in the representation of the instantaneous eigenvectors 
(33) reads: 

_l 



U{t,t') = 



(cos(e(t - t')) - ir]sm{e{t - t')))e^/'^'^^*-''^ -i/?sin(e(i - i'))e^/2¥'(t+*') 

-i(3sm{e{t - t'))e-V2¥'(W) (cos(e(t - t')) + iTjsm{e{t - t')))e-'/2'^(*-*' 



(39) 



where e = (^ + ^ cos'!?)^ -I- and s = ^siniS. 
Also: 

Wo r+^cosz? 
2e e 

This is the propagator in the adiabatic basis. Therefore, 



in order to switch to the fixed spin basis, one should write 
Usptnit,t') = B{t)U{t,t')B\t'), where B{t) is given by 
Eq. . The four elementary stretches imply the follow- 
ing substitutions in the propagator of Eg. ([39)) : 

u^) forward orbiting in the upper arm of the ring : 

ip{t) — LOot and i9 — > 
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u^) backward orbiting in the upper arm of the ring : 
ip{t) = TT — ujot and 1} — *■ 

d^) forward orbiting in the lower arm of the ring : 
ip{t) = 2tt — uJot and 1} — > — i?. 

backward orbiting in the lower arm of the ring : 
ip{t) = TT + UJot and 1!) ^ 



VII. THE CONDUCTANCE 

In this section, we derive the DC conductance Q across 
the ring, at the Fermi energy. Within Landauer's ap- 
proach, Q is given by 



e 

'h 



Y,\A{a;a'\Eo)f 



(40) 



We will here consider the dependence on the external 
magnetic field {(t>/4>o) and on the spin-orbit strength 
ksoR^^ both in absence and in presence of dephasing 
at the contacts. The various amplitudes in Eq. (j40ll have 
been numerically computed, as discussed in Sec.s(III-VI). 
In Fig. (HI we show the magnetoconductance of the ring in 




FIG. 4: (color online) a) Magnetoconductance of an ideal ring 
as a function of the magnetic field 4>/4>o- The Fourier trans- 
form (FT) of the magnetoconductance (in the inset) shows 
only the AB peak at freq. 4>o^ ■ b) Magnetoconductance of 
a realistic ring in which the S matrix of Eq.([71) regulates the 
scattering of the electron at the leads. As a consequence of 
the backscattering, the FT of the magnetoconductance (inset) 
shows higher order frequencies. 

the absence of RSOI [ksoR = 0). In panel a) of Fig.®, 
only the path of Fig. ^ la) has been considered, i.e., full 
transmission across the ring is assumed, as it would be 
the case for ideal coupling to the leads. The correspond- 
ing Fourier spectrum is showed in the corresponding in- 
set. The well known Ahronov-Bohm sinusoidal pattern 
implies that just the fundamental frequency 0q ^ appears. 

To make the model more realistic, we allow for higher 
order looping of the electron within the ring. In Ref.^^, 



only the paths of the kind of Fig.([l]2a)[l]2b)) were in- 
cluded. Here, we consider also the paths of the kind 
of Fig.(l2|) in which the electron can be backscattered at 
the leads. We use here f = in the scattering matrix 
between the arms and the leads, which means that no 
backreflection in the incoming lead is present. The mag- 
netoconductance of the ring, in this regime, is showed 
in Fig. (IDd). Because of the inclusion of time reversed 
paths (TRP) within the ring, we see that higher order 
frequencies appear in the Fourier spectrum. In particu- 
lar, the inset shows a peak at 2/(j)Q which is the signature 
of weak localization™. We also include some dephasing 
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FIG. 5: (color online) Conductance as a function of 4'/4'o (left 
panels) and ksoR (right panels) for ideal (black curves) and 
realistic (red curves) contacts. An increasing amount of de- 
phasing at the contacts is also included: from top to bottom: 
^ = tt/3, n, 2-K 



due to diffusiveness in the contacts by adding a random 
phase z e (— C, C) for each scattering at the leads. In 
Fig-®, we report the conductance G vs. (j)/4>o , with 
ksoR = (left panels) or vs. ksoR with 0/00 = (right 
panels). These are averaged over N = 1000 realizations 
of disorder, and plotted increasing the window of phase 
randomness (^ — 7r/3,7r, 27r from top to bottom). The 
black curves refer to the case of Fig(|3] a)) (ideal con- 
tacts) while the red curves refer to the case of realistic 
contacts (Fig^b)), with f = 0. 

By comparing the top left panel of Fig. ([5]) with Fig. ([4]), 
we see that the ring is rather insensitive to small de- 
phasing at the contacts. By increasing the amount of 
dephasing (middle and bottom left panels in Fig. |5P ) we 
find that the sensitiveness is larger in the case of real- 
istic contacts. This is due to the fact that for realistic 
coupling, the electrons in the ring can experience higher 
order paths, since it scatters with the leads many times. 

In the right panel of Fig. © , we plot the DC conduc- 
tance vs. ksoR at 4>/4>o = for both ideal contacts 
and realistic contacts (and f = 0) (black and red lines 
in each box), with an increasing phase randomization 
(boxes from top to bottom with ^ = 7r/3, tt, 2tt), averaged 
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over TV = 1000 realizations. In the case of ideal contacts 
and little dephasing (top right panel black curve), we 
see again the quasiperiodic oscillation of the conductance 
reproducing the localization conditions at the expected 
values of fcso- B^^'^^'^'^'^^'^'' - When including higher or- 
der processes, interference effects give rise to a slightly 
different pattern. In the case of realistic contacts, we 
note that the device is seriously affected by dephasing, 
mainly because including the TRPs contributing to the 
transmission amplitude increases the number of scatter- 
ing processes at the leads. Indeed, when the dephasing 
is quite large, it gives rise to random oscillations that are 
not averaged out, so that they wash out the conductance 
oscillations. The effect takes place for C ~ tt when TRPs 
are included, in contrast to ^ ^ 27r when the TRPs are 
absent. As regular magnetoconductance oscillations are 
experimenally observe d^^'^^i^^'^"^i^^ with little precentage 
of contrast between maxima and minima, we conclude 
that, in real samples, dephasing is ubiquitous. 



VIII. SPIN TRANSMISSION 



and the final transmission of the spin will result to be 
partly out of the x — y plane. On the contrary, when 
4'/4'o ^ ksoR the competition of the Zeeman and the 
RSOI induces a coherent rotation of the spin while the 
electron travels along the ring. Fig.® shows that, when 
4>l 4'o ^ ksoR, the spin is moved significantly out of the 
x — y plane, consistently affecting the AB oscillations. To 







H 


\r^} 


w 








„R=4: 



In this Section we calculate the rotation of the spin 
of the electron transmitted through the ring. We first 
consider an incoming electron beam with in-plane spin 
polarization (let's say, polarized along the x direction ). 
The spin rotation is measured by calculating the average 
value of the outgoing spin: 



I 



(41) 



Since in the previous Section we have shown that higher 
order looping just adds subleading higher order harmon- 
ics to the conductance, here we focus on the case of 
ideal contacts, that is, we include in the calculation only 
paths as the ones of Fig.lUla). The in plane polarization 
can be considered as a superposition of equal weighted 
z-polarized spin components. In the absence of RSOI 
{ksoR — 0), opposite spin polarizations do not interfere 
with each other. As a consequence, the total expected 
{Sz) component keeps zero at the exit. Fig.® shows the 
magnetoconductance for increasing ksoR, and the cor- 
responding expected spin component polarized along the 
z— axis at the exit of the ring. The Zeeman term is on 
the diagonal of the spin Hamiltonian of Eq. (pij) and op- 
erates to keep the z— components of the spin polarization 
fixed, while the RSOI is offdiagonal and tends to favor 
inplane spin polarization. This implies that when the 
magnetic field increases {4>/4'o >> ksoR) the spin po- 
larization gets frozen to the incoming polarization. This 
is the reason why,in Fig.®, at high fields, the trasmit- 
ted spin polarization is inplane. Incidentally we observe 
that this result should not be expected in real systems in 
which spin relaxation can occur due to electron-phonon 
interaction or other mechanisms as hyperfine interaction 
with nuclear spins. Spin relaxation would induce flipping 
of that spin component that is energetically unfavourable 



FIG. 6: Magnetoconductance and expectation value of the 
outgoing z spin component for an incoming spin in the x 
direction at different values of the RSOI strength {ksoR = 
0, 1, 2, 4 indicated in the pictures). 



better understand what happens when 4>/(j)o ^ ksoR, 
we isolate the spin "up" polarization for the incoming 
particle in what follows ((*j„| l^-m) / {'i!^n \^^n) = 1) 
and we separately plot in Fig. ([7]) the two contributions to 



the conductance G 



up— up 



[full line) and G, 



down— up 



{dotted 

line), for opposite polarizations of the outgoing electron. 
Gup-up is the contribution to the conductance arising 
from the particle flux that mantains the same polariza- 
tion at the exit as the incoming one, while G down— up 
refers to a particle flux having the opposite polariza- 
tions at the exit with respect to the one at the en- 
trance. Of course, when ksoR — 0; the electron spin 
is in the "up" direction for any 4'/(l)o. When both RSOI 
and magnetic fleld are present, with 4'/4>o » ksoRj 
the spin polarization is steadly in the z— direction, ex- 
cept for sharp reversals at flux values (f>om/2 (m inte- 
ger). However, G down-up is vanishingly small at these 
places, together with G 



up— up 



Therefore the conductance 
vanishes at these points anyhow, and the transmitted 
spin polarization is fully up, except for these points. On 
the contrary, in the parameter intervals characterized by 



ksoR, both Gil 



and G, 



down— up 



are non van- 



ishing (see Fig[7l), so that the spin is rotated at the exit, 
with nonvanishing transmission amplitude. 

We now examine in more detail the expected depen- 
dence of the outgoing spin polarization on the RSOI, at 
zero magnetic field. For an incoming electron polarized 
with spin "up", Fig.([5]Ja)) shows that a large enough 
RSOI produces a rotation of the spin which points down 
at the exit for any ksoR value. Meanwhile, the total con- 
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FIG. 7: Separate contributions to the magnetoconductance 
for different outcoming spin polarizations, Gup-up {full line) 
and Gdown-up {dotted line ) compared to the expected value 
of the outgoing z spin component. The incoming spin is polar- 
ized "up" . Different values of the RSOI strength are reported 
{ksoR = 0, 1, 2, 4 indicated in the picture ). 



ductance oscillates with ksoR- This finding also appears 
in Refi^ and is quite remarkable, because it is the con- 
sequence of the interference between the two arms of the 
ring. In order to point out the role of quantum interfer- 
ence, we discuss here, for comparison, what happens by 
transporting a spinful electron along a single arm of the 
ring (the upper one). When just one arm is considered, 
(see Fig(|5jD)) spin polarization oscillates, as a function of 
ksoR^'^ 1 while the conductance is always unitary, because 
of the conservation of the particle flux. In the limit of 
large ksoR the propagation amplitude for the travelling 
electron acquires a simple analytical form. From Eq. (j39p 
we see that, in the representation of the instantaneous 
spin eigenvectors, the spin propagator at the exit time i/ 
(with ojotf — T^), in this limit is diagonal: 



(i/,0) = 



ie-'''^ 
-ie"'' 



(42) 



so that the spin appears not to be rotated at the exit 
in the rotating reference frame. However, in order to 
move from the representation of the instantaneous spin 
eigenvectors to the reference basis, one has to perform 
the transformation with the unitary B matrix of Eq. (IM)) 
, with I? = 7r/2. The spin part of the propagator for 
an electron travelling into the upper arm (upper path of 
Fig.dUla)) is then: 



i cos(7r7) sin(7r7) 
— sin(7r7) —i cos(7r7) 



(43) 



If we inject up electrons in the upper arm only, the expec 
tation value of the the outcoming Sz defined in Eq.lHT] 



is: {Sz} = cos(27r7) = cos{nksoR) (note that for large 
enough 7 this result well approximates the red-full line in 
Fig.lIHjD)). The conductance is unitary, G — 2e^//i (the 
black-dashed line in Fig.([S)3), as no interference takes 
place. 

We now go back to the transmission along both arms 
simultaneously and examine the resulting interference. 
According to the rules given after Eq. ipS)) . in the same 
limiting case as above, the propagator accounting for 
transmission of incoming up spins across the ring is: 



spin 



2sin(7r7) 
— 2sin(7r7) 



(44) 



so that the spin at the exit is reversed. In fact, in the 
expectation value of Eq. (|¥T|l . the oscillations in the nu- 
merator compensate those in the denominator, eventu- 
ally giving {Sz) = —1 (for large enough 7 this result 
well approximate the red- full line in Fig. ((5^)). The con- 
ductance however oscillates according to G/{2e'^/h) = 
2 sin^(7r/c5oi?/2), as plotted in the black-dashed line in 
(FiglSi)). ^ 

It is quite remarkable that this result is only found 
at zero magnetic field. Indeed, no matter how small B 
is, the time reversal symmetry is broken and the spin 
oscillates with ksoR (see Fig.(j9la). However, for very 
small magnetic field these oscillations are confined close 
to special values of ksoR — 2^(Z integer) and display a 
Lorentzian shape around these points. The role of the 
magnetic field is to make the oscillations broader. 

To summarize, there are two limiting conditions in the 
outgoing spin polarization, for incoming "up" spin po- 
larization: a) the zero RSOI which leaves the incom- 
ing spin polarization unchanged, provided no relaxation 
takes place; b) the zero magnetic Eux case in which the 
RSOI produces a flip of the spin at the exit. It is in- 
teresting that when the flux is an integer number of 
flux quanta the crossover between case a) and case 
b), with ksoR increasing from the value zero to values 
ksoR >> 4'l4>o is rather sharp. This is shown in Fig. 
(Ob) where the expectation value of the outcoming spin 
is plotted vs. ksoR for different integer values of 4'/4'o- 
For values of 4i/(j)o » ksoR the outgoing spin polar- 
ization is the same as that at the entrance, ( up in the 
picture). On the contrary, by increasing ksoR £^t non 
zero B field, we see again a pattern similar to the one of 
Fig. (IHa), but shifted to higher values of ksoR- 



IX. CONCLUSIONS 

To conclude, we have employed a path integral real 
time approach to compute the DC conductance of a 
ballistic one dimensional mesoscopic ring in both exter- 
nal electrical and magnetic fields orthogonal to the ring 
plane. The spinful electron experiences a Rashba spin or- 
bit interaction and the Zeeman term. We employ a piece- 
wise saddle point approximation for the orbital motion, 
but we implement the full scattering matrix at the leads 
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FIG. 8: (color online) Conductance {broken line) and spin 
polarization {full line) of the outgoing electron for B = 0, as 
a function of the RSOI. The spin of the incoming electrons 
is polarized "up" : a) The ring case with ideal contacts, b) 
a single wire of the same length and curvature as one of the 
ring arms. 















0-^) (j) 






FIG. 9: (color online) Expectation value of the outgoing spin 
for incoming spin "wp" polarized electrons as a function of the 
RSOI for different values of the magnetic flux : a) <j)/4'o zero 
or very small; b) increasing integer values of (f)/(j>o. 



and sum over all the possible higher order paths up to 
convergency of the result. Our approach goes beyond 
other recent semiclassical calculations. Our theory is 
nonperturbative and separates the adiabatic spin dynam- 
ics from the non adiabatic one by using the rotating frame 
for the spin travelling around the ring. In practice, we 
diagonalize the time dependent spin Hamiltonian in the 
representation of the spin eigenvectors of the instanta- 
neous Hamiltonian. This allows us to explore a wide 
range of Hamiltonian parameters, ranging from the limit 
of strong magnetic field and weak Rashba SOI to the 
opposite case. In both extreme regimes our piecewise 
saddle point approximation is very efficient as quantum 
fluctuations with flipping of the spin has little influence 



on the orbital motion. This is also seen from the num- 
ber of paths required to gain full convergency. As ex- 
plained in Sec. V the separation of adiabatic from non 
adiabatic spin dynamics shows that in the intermediate 
regimes our approximation is less justiflcd, but neverthe- 
less, the results it produces seem to be in rather good 
agreement with recent numerical calculations2Si22ii2ii2^i22 
and experimentsiiii^ii^iii^ii^. When we include also time 
reversed paths, the Fourier transform of the magneto- 
conductance shows the typical (f)o/2 peak due to weak 
localizationi^.This would be the only surviving contri- 
bution if an ensemble average or different rings were 
measured?^ . 

We have also allowed for nonideal couplings between 
ring and leads as we account for dephasing effects due to 
diffusiveness at the contacts. The results satisfactorily 
compare with experiments where the contrast between 
maxima and minima in the interference fringes is always 
few tens of percentage of the background DC signal. 



Acknowledgments 

We acknowledge financial support by the Italian Min- 
istry of Education (PRIN) and by the CNR within ESF 
Eurocores Programme FoNE (contract N. ERAS-CT- 
2003-980409). 



APPENDIX A: MOTION OF A CLASSICAL 
SPIN IN A ROTATING MAGNETIC FIELD 

In this appendix we derive the classical equations of 
motion for the spin from the Lagrangian in Ea. p4p . by 
assuming for the orbital coordinate the saddle point so- 
lution ip = constant. Once the orbital motion is dealt 
with in this way, the Lagrangian for the spin degrees of 
freedom is given by (besides a constant contribution) 



£[6,$,$]/^ = 



1 — cos 9 



(Al) 



where the effective time dependent magnetic field is 
B = {B+,B^,B^) ^ (7(^e*'^,7(y9e-*'^,-Wc/2). To derive 
the equations of motion from a variational principle, we 
write the Berry phase term in the total spin action as 



Se 



*(*/) 



dt 



*(0) 



1 — cos O 



d$ = / sin e rfe A d$ 



(A2) 

where S is the spherical triangle with vertices given 
by the north pole on the sphere and by the points with 
coordinates (9(0), *(0)) , (9(i/), $(t/)). Let {t,u) be 
a parametrization of the spherical triangle, such that 
5^(^,1) = S{t), and S{t,{)) = (0,0,1). Thus, one may 
rewrite the action Sb in Ea. (|A2[) as 



Sb - 



.T .1 


'dS 


dS' 


/ dt 1 duS ■ 






10 Jo 


9i 


du_ 



(A3) 
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To derive the equations of motion, we consider a vari- 
ation S{t,u) S(t,u) + SS{t,u) such that S{T,u) and 
S{0,u) are "locked", that is, SS{0,u) = dS{T,u) = 0. 
Since [S{t, u)]^ = 1 Vt, u, one gets s'- ff = ^- f| = 0, As 

a consequence, ^ x ff is paraUel to S. As SS ■ S ^ 0, 
this imphes that 



dt / du6S 
Jo 



dS_ dS_ 

dt du 







(A4) 



Thus, by integrating by parts we see that the only 
nonzero variation of Sb is given by the boundary term 



5Se 



dt SS{t) 



Sit) X 



dS{t) 
dt 



(A5) 



where we have used the fact that S{t,l) — S{t). By 
equating to zero the total variation of the action, one 
obtains 



dt 



(A6) 



that is, the classical equations of motion we used in sec- 
tion IV. To show that Eqs. (jA6p . when the spin compo- 
nents are expressed in polar coordinates, are equivalent 
to Eas. (|17ll8p . let us set uj = ^ip and = — cjc/2. Also, 
we define 



S+ — Sx 



iSy 



Sx iSy 



B+=Bx + iBy 

B-^Bx~lBy . 



In terms of the new variables, the equations of motion 
are given by 



dS^_ 
dt 

dS: 



-{B+S- -B-S^ 



dt 



i{B,S, 



B+S, 



dS_ 

~dr 



-i{B,S--B-S,) 



(A7) 



dm{t) 
dt 



i (yi ~ ip) p{t) ~ 2 iuj Sz 



dt 



dSz iuj , . , . „, 

^ = -y-W (A8) 



where 



p{t) = S+e^-'f''^ + S^e^'f'^ 
m{t) = S+e^-"^^ - S-e'-"^^ , 

and 1 = 4|S'p = 45^^ + (p2 - m^). By introducing b = 
CI ~ If, we obtain: 



d(m(t) +pit)) 
dt 
dSzjt) 
dt 



ib{m{t) +p{t)) - 2iujSz{t) (A9) 



-iuj m{t) 



(AlO) 



Resorting to the polar coordinates (8, $) for the spin S, 
we get: 



e cos e + i $ - sin e 



icosQaifi = O(All) 



e 



uj sm X 



sine = .(A12) 



Ea. (|XT2)) is the same as Ea. (|17p . The real part of 
Eq. ljAlip is proportional to the imaginary part: both 
give Eq.((T8]) when equated to zero, which completes the 
proof. 



APPENDIX B: THE SPIN PROPAGATOR 

In order to find the propagator of the Berry Hamil- 
tonian Ha of Eg. ([57)1 . we solve the system of differential 
Eq. ([551) I ill the representation of the instantaneous eigen- 
vectors: 



I — 

dt 



r + ujo sin^ ^ 
ujo sin ^ cos |-e*'^ 



,sin|cos|e-»"°* 

• 2 i9 

— r — ojo sm ■2 



U4 
7i_ 



(Bl) 



To solve Eq. (|Bl|) . first of all, we switch to a time- By setting 
independent coefficient matrix by defining: 



we define the matrix T through 

y+ \ _ f e+'^* \ f u ^ 



Y ^ 


I- 


; W = 













y- J \ e-'^* J J ■ ^^"^^ Y = TW, W = T-^Y. (B3) 
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Eqs. ljBip now read: 



which, by defining V = P ^ Y, becomes 



= ir- ^^cosmy+it) + ^ sin^y^it) 



^^it) = +Y sindy+it) + {-r + ^cosmy^{t) {B4) 



Now we define r' — r — ■ 
in matrix form we have: 



Y cos-i? and s — ^ shi'd, so that 



r s 
s — r' 



y+ 



(B5) 



dt \ y 

in a compact form we can rewrite the last equation as: 
. d 



(B6) 



which defines the matrix C. 

We now decouple the previous system of equation by di- 
agonalizing the matrix C. Its eigenvalues are A = ±e = 
±\^r''^ + and the matrix that diagonalizes C is 



dt V -e ' 



or, in matrix form 



V{t) = S{t~t')V{t') 



(BIO) 



Its formal solution is: 



(B12) 



Now we apply inverse transformations, in order to ob- 
tain the full Schrodinger propagator, that is the matrix 
transformation between {W{t') and W{t)). 



P 



1 

e— r 



r — £ 
s 

1 



(B7) 



W{t) = T-^{t) P S P-^ T{t')W{t') 



Its inverse is 



2£(£-r') 2e(e-r') 

__(C-V)£_ 

2e(£-r') 2e(e-r') 



Eq. (jB6[) now reads 
. d 



i—P-'Y ^ P-'CPP-'Y 
dt 



(B8) 



(B9) 



where P is time independent. The full evolution operator 
in the adiabatic basis is: 



U{t,t') = T-^{t) P S p-^T{t') 



(B13) 



By performing all the matrix products, we obtain the 
result of Eq. ((39|) . as given in the text. 
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